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ABSTRACT 


The asymptotic relation between the pressure with exchange in the 
Thomas -Fermi-Dirac atom model and the corresponding pressure 
without exchange on the Thomas -Fermi model for the case of low atomic 
number considered in the first paper of this series is generalized by 
including a correction corresponding to the first-order Coulomb con- 
tribution to the energy of the atom, rather than the kinetic energy of 
the electrons alone. All terms yielded by the asymptotic solution of 
March for the Thomas-Fermi-Dirac equation in the limit of small atom 
radius are included, and the extension of results on the basis of the 
general asymptotic solution of Rijnierse in this limit is indicated. The 
expression for the radius of an isolated Thomas-Fermi-Dirac atom 
for low atomic number obtained previously is extended to include a 
term corresponding to the Coulomb correction of first order, and com- 
pared with a similar result of Jensen. Asymptotic relations between 
thermodynamic functions (the energy, the compressibility, and param- 
eters associated with the energy and equation of state) on the two models 
for low atomic number as given previously are extended similarly. 
Comparison of the extended asymptotic relations between the pressures 
on the two models is made with numerical results obtained by direct 
solution of the Thomas-Fermi-Dirac equation. 
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I. INTRODUCTION 


An asymptotic relation expressing the pressure with exchange on 
on the Thomas -Fermi-Dirac atom model as a function of the corre- 
sponding pressure withojat exchange on the Thomas-Fermi model was 
obtained by Gilvarry^ on the basis of the method of March^ for solving 
the corresponding differential equations of the models for an infinitesi- 
mal atom. In a preceeding paper by Gilvarry, Hartle, and March, ^ 
hereafter referred to as I, it was shown that the result of Gilvarry 
follows by a direct physical argument, without the necessity of recourse 
to solution of differential equations. This possibility arises because the 
relation between the pressures correspnds to the limiting situation at 
small volume, where only the kinetic energy of the electrons distributed 
uniformly throughout the atom need be taken into account (whether for 
high pressure or for vanishing atomic number). 

It was pointed out in I that the physical method presented there 
could no be used to generalize the results of Gilvarry. The purpose of 
this paper is to obtain the generalization of these results by taking into 
account all the terms in the asymptotic solution of the Thomas -Fermi- 
Dirac equation for small atom radius which can be determined by the 
method of March. Physically, this procedure corresponds to consid- 
ering in the relation between the pressures with and without exchange 
the effect of the Coulomb interactions to first order as well as the 
effect of the kinetic energy of the electrons. In addition, the expression 
obtained in I for the radius of a Thomas-Fermi-Dirac atom in the limit 
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of vanishing atomic number will be correspondingly generalized, and 
the connection with prior results of Jensen^ will be discussed. Results 
for thermodynamic parameters as generalized in the case of low atomic 
number also will be eonsidered. 

Only a restricted number of terms in the solution of the Thomas- 
Fermi-Dirac equation for the limit of vanishing atom radius can be 
obtained by the method of March, because of the appearance of diver- 
gences. In particular, only terms through quadratic in the exchange 
parameter can be found. However, Rijnierse® has developed a method 
of obtaining systematically any desired number of terms (and of arbi- 
trary order in the exchange parameter) in the solution of the Thomas - 
Fermi-Dirac equation for the case of small atom radius. The extension 
of results on the basis of the general solution by Rijnierse will be 
indicated. 

The treatment is restricted to the case of zero absolute tempera- 
ture throughout. The atom is considered spherically symmetric. 

II. ASYMPTOTIC RELATIONS 

The Thomas -Fermi-Dirac equation® is 

d^'/'/dx^ = x[('///x)^^^ + e]^, (1) 

where e is the exchange parameter ! 2>Z'^ ) Z for an atom of 

atomic number Z, and x is a dimensionless variable connected with the 
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ridial distance r in the atom by r = /xx , in which yu. is the scale length 
( 97r2/i 28 Z ) 1/3 in atomic units. The Thomas-Fermi-Dirac function 
i// (x) for a neutral atom satisfies the initial and boundary conditions 

«0)=1, [dy,/dx]^^^^ (2a, b) 

respectively, where is the boundary value of <// corresponding to the 
radius Xj^ of the atom- When e vanishes, Eq. (1) reduces to l3ae Thomas - 
Fermi equation d^^/dx^ = 03/2/^ 1/2 ^ where the Thomas -Fermi 

function 0 is subject to initial and boundary conditions exactly analogous 
to those of Eqs. (2). 

Asymptotic solutions for the Thomas-Fermi-Dirac function in the 
case of small atom radius have been obtained by March and Gilvarry 
by means of a method introduced by the former. The procedure 
consists in expanding 0 in a Taylor series about the atom boundary, 
imposing the initial and boundary conditions, and obtaining systemati- 
cally the partial sums 0^ of degree n = 1 , 3 , 5 , * * * in x - Xj^ . The 
first three partial sums are given recursively by 

0j = 33/3 (x/x^) _ 2 • 3i/3(x/xj^) e + xe^ (3a) 

03 = 0j + 1 - £ (x/xj^) +i (x/xj^)3 

'/'s = '/'s - (3/10) x/x^. 


(3b) 

(3c) 
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These results are valid in the limit of vanishing atom radius . 

It was shown by Gilvarry and in I that the Thomas -Fermi -Dirac 
function corresponding to the limit of a uniform distribution of electrons 
confined to the atom volume is the partial sum 0^ of Eq. ( 3 a) and must 
include terms linear and quadratic in e . Of the three terms, the initial 
one arising from the chemical potential of an electron is the dominant 
one and the last two correspond to the effect of exchange. The first 
term dependent on e was obtained by March and the second by Gilvarry. 
The original results of Gilvarry and those of I follow from the boundary 
value of , upon which the pressure depends. This partial sum 

corresponds to the effect of the kinetic energy of the electrons only. 

The second partial sum 0 ^ of Eq. ( 3 b) differs from 0^ by the last 
three terms, which represent the effect of the variation of the potential 
in the atom computed on the assumption that the electron density is 
uniform. These terms were obtained first by Gilvarry, by means of a 
physical argument. However, one notes that the three terms cancel 
for X = Xj^ , and hence the boundary values 0^ and 03 corresponding 
to the first two partial sums are equal. It follows that the associated 
pressures agree and hence the effect on the Thomas-Fermi-Dirac 
function of the variation of the potential implied by a uniform electron 
distribution does not change the results of I inferred by considering 
the effect of the kinetic energy alone. This fact was noted by Gilvarry. 
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The partial sum 1//5 differs from 4^3 by a term — (3/10) x/xj^ which 
arises from the contribution of the electrostatic electron-'nucleus and 
electron-electron interactions in addition to the kinetic energy, producing 
nonxiniformity of the electron distribution in the atom. Since this 
correction to does not vanish at the atom boundary, the corresponding 
boundary value 1 //^ does not equal and the associatedipressures 

differ. The_ob.ject of this paper is to obtain the modificationdo the 
results of I aiasing from inclusion in the Thomas -Fermi-Dirac function 
of the effect of the term in question, arising from the electrostatic 
interaction when the electron distribution is not uniform. This term 
representing the Coulomb correction on the boundary vaihe of the 
Thomas-Fermis- Dirac function was obtained first by March. 

Note that S<q. (3c) for contains all terms in^p which can be de- 
termined by the method of March. Terms of higher order are given by 
Rijnierse. 

In what follbws, it is necessary to consider only the boundary 
value 

'/'(j = - 2*3^^^ e + Xj^€^ -3/10 , 

corresponding to \p^, since the pressure depends only upon. boundary 
values. In this expression, the first three terms correspond to the re- 
sults of Gilvarry and in I; the last term arises from the Cbulomb cor- 
rection and yields the extension of the results of I obtained’in this paper. 
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The pressure with exchange corresponding to a Thoma s- Fermi - 
Dirac atom at zero temperature is given by 



where e is the electronic charge. When e vanishes, this expression 
reduces to the corresponding pressure P of a Thomas -Fermi atom 
given by 


in terms of cf:)^ , the value of cp at the atom boundary. Upon substituting 
Eq. (4) into Eq. (5) and retaining terms of significant order, one obtains 

Pex = Po (7) 

where Pg is the Thomas-Fermi pressure in the limit given by® 

Pg = (2V157 t 3) (2/7^2)173 (eVa^) x"® (8) 

in terms of the radius Sg of the first Bohr orbit for hydrogen. 

The expression for the pressure p^^ of Eq. (7) was obtained by 
March2 from a form for \jj^ corresponding to Eq. (4) with the exception 
that the quadratic term in e was ommited; Gilvarry pointed out that 
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March's result for is valid through terms in e higher than were 
included explicitly. The leading term of Eq. (7) is just the contri- 
bution to the pressure from the kinetic energy of a uniform distribution 
of Z electrons in the atom volume. The second term is simply the 
exchange correction to such a gas of free electrons. The laS^t term can 
be identified as’ the sum of the pressure contributions p^^j and p^.^ from 
the electron-nttcleus and electron-electron interactions, respectively, 
calculated on the assumption of a uniform distribution of electrons. 

It follows from Eq. (7) that the Thomas -Fermi-Dirac pressure 
p^^ can be expressed as a function of the limiting form Pq of Eq. (8) 
for the Thomas-Fermi pressure by the asymptotic relation 


p^^ = Po [ 1 - (1+32/3 e-VS)(P/Po)i/5] (9) 

in terms of the constant pressure (independent of Z ) given by 

Pg = [5^/(3-2^°w7)]e2/a4, (10) 

as was noted by Gilvarry. The fact that the Thomas -Fermi-Dirac 
pressure p^^ of Eq. (9) includes the effects to first order of the electron- 
nucleus and electron-electron interactions raises the possibility that 
Pg^ can be expressed as a function of the corresponding Thomas -Fermi 
pressure which includes the effects of these interactions. One obtains 
the latter pressure including such interactions to first order by equating 
e to zero in Eq. (7) to yield 
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Pm = Po , 


(M) 


where is the pressure derived by March. ^ Factoring (l-3^^^Xj^/4) 
from the expression in square brackets of Eq. (7), one obtains 

P«x = Pm tl- (5ex^/4-3i/3) + , (12) 

which differs from Eq. (7) in that is replaced by Pjj , which contains 
the Coulomb correction. 

The terms neglected in the square brackets of Eq. (7) are of form 
0( ). The perturbation analysis of March leading to this equation 

cannot yield the coefficient of the termO( ); however, this coeffi- 
cient can be determined by consideration of the general form of the 
Thomas -Fermi-Dirac energy . The energy to which Eq. (7) corre- 
sponds is^ 

’Jex = ""o [l-(5£x/2-3i/3)_ 

where Uq is the Thomas-Fermi energy associated with p^ . One can 
infer from this expression and the results of Rijnierse that the general 
asymptotic form of u^^ for the case of small x^^ must be 

00 

u = u 0 y b X? ( ^ “^) 

ex 0 / n b 

n = 0 
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witli the coefficietits dependent on e in general. Since Uq varies as 
x“ , the corresponding pressure is 


Pex = PO 


L 

n = 0 



( 15 ) 


and one notes the complete absence of a term of form Pg 0( Xj^^ ), since 
the corresponding coefficient vanishes when n = 2 . This conclusion is 
consonant with the general solution of Rijnierse for p given by Eq. 

(Al) in the Appendix. 

One notes that the factoring operation leading to Eq. (12) introduced 
a term of order into the expression in square brackets. It follows 

from the preceding argument that the corresponding coefficient is exact, 
since the form of p^^ preceding the factoring operation contained no 
such term. This fact is important in what follows, since it will be 
demanded that the final asymptotic relation for p^^^ vanish when the 
radius of the atom has the value ^ corresponding to the isolated 
Thomas -Fermi-Dirac atom. Terms of order must be retained 

in the asymptotic expression for p^^^ , since they are of the same order 
as Xj^ 2 when Xj^ has the value ^ corresponding to the isolated atom, 
as will appear. 

Equation (12) can be put in the form 

Pex = Pm fl- (5eXj,/4-3i/3) (1 + ^^ 20 )- £x2/4] (16) 
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by adding and subtracting ex^^/16» the limit of small Xj^ at a 
fixed value of Z, the Thomas-Fermi pressure p approaches asymptotically 
the values of P^, and Pq . However, for small , in general, Pj, should be 
a better approximation to p than is . Hence, for the limit in question, 

Eq. (16) can be written with suppression of the subscript M onp as 

P„ = P [1- (P„/p>‘'* (Pi/p)2«l (17) 

by substitution for e and elimination of x^^ with use of Eq. (11). The 
constant pressure is given by 

Pj = (4/5)5 ( 3 ^ 2 / 2)576 ^ (18) 

in terms of Pq of Eq. (10). 

Equation (17) represents the proper asymptotic relation between 
the Thomas-Fermi-Dirac pressure p^^^ and the Thomas-Fermi pressure 
p when the^ effects on the pressures of not only the kinetic energy of the 
electrons but also their Coulomb interaction to first order are taken 
into account. On the same basis as in the argument of Gilvarry and as 
given in I, the result holds in the first instance for infinitesimal volume 
and should be valid in the two limits corresponding to infinite pressure 
for arbitrary atomic number and to vanishing atomic number for arbi- 
trary pressure. In general, Eq. (17) should hold for a range of atomic 
number at low Z, rather than simply in the asymptotic limit Z-‘0. 
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One notes from Eq. (17) that one obtains Pgx = P in the limit p-*oo 
at constant Z, where the contributions from the exchange and Coulomb 
energies vanish. In contrast to the limit at fixed Z corresponding 

to consider the limit Z-* Oat a fixed value of the Thomas- Fermi 

pressure p . Forthis’case, the asymptotic expression of Eq. (8) implies 
Xj^-*0 as Z-‘0, independently of pressure P. The last terra in Eq. (17) 
representing a contribution from the Coulomb energy vanishes for 
the limit in question and the resulting expression 


Pex = P (19) 

is that obtained by Gilvarry and in I. 

Further generalization of the asymptotic relation of p^^^ to p is 
carried out in the Appendix by using the general expression for P^^j 
obtained by Rijnierse. In this instance, for sufficiently small , Eq. 
(A5) can be written as 


p^^ = p [l-(P/p)l/5_Z2/3 (P/p)2/5_ Z4/3 (P^/p)3/5] (20) 


in terms of the constant pressure 


p2 = [(2/5)'* 


( 21 ) 



This procedure yields a result that should be valid over wider ranges 
of pressure than the corresponding Eq. (17). One notes, of course, that 
the same asymptotic results follow from either Eq. (17) or Eq. (20) 
when one takes the limits P-<oo at fixed Z or Z-*0 at fixed P . 

A maximum atom radious 2 value of Z corresponds to 

the isolated Thomas -Fermi-Dirac atom with zero pressure on its 
boundary. The boundary value '/'j, ^ of '/'j, in the case p^^^ = 0 must 
satisfy the Jensen boundary condition 

'/'b,z'’'b.z = (22) 

which constrains 2 1° ^ function of Z . The Jensen condition 

conjoined witli ijj^ of Eq. (4) yields a quadratic equation for ^ with 
a corresponding ambiguous sign in the root. Eq. (17) requires P>I®q 
and since p( ) is a monotone-decreasing function,^ the minumum 
root must be taken. Hence, the radius ^ of the Thomas -Fermi-Dirac 
atom on the basis of Eq. (4) is 

32/3 g- A ( 23 ) 


b,2 


.31/3 ^-1 (1. 


or, equivalently. 




1-- (96 7t2)1/3z2/3 


(24) 
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where terms of significant order have been retained. One verifies 
directly that 2 2 same order, as stated 

previously. 

Direct substitution of the result of Eq. (24) into of Eq. (17) 
yields 

Pex(^b.z)=0 

to the order of terms considered. Thus, Eq. (24) represents the radius 
of an isolated Thomas -Fermi -Dirac atom as a function of atomic number 
in the limit Z-»0 , in conformity with the assumption that Xj^ be small. 

The leading tearm of Eq. (24) , obtained in I, corresponds to the effect 
produced by the exchange force which tends to contract the electron 
cloud (through its effect on the electrostatic potential energy). The 
second term corresponds to a contribution from the Coulomb interaction, 
in conformity with the fact that the total electrostatic attraction of the 
nucleus on the electrons is greater than the repulsion of the electrons 
on themselves. To the isolated Thomas -Fermi -Dirac atom with radius 
Xfa z given by Eq. (24) and meeting Eq. (25), a compressed Thomas - 
Fermi atom of the same radius exists for which the pressure P satisfies 


pK.z) =1^0 [l+(16/5) (37rV2)i/3 Z^/^] 


( 26 ) 
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This expression generalizes a result of I to low atomic number, since 
it yields the physical interpretation of 1®^ given by p(xj^ 2 ) ~ ^0 
limit Z -» 0 . Notice that as Z increases, the pressure at the boundary 
of the isolated Thomas -Fermi-Dirac atom remains zero, but that for 
the associated Thomas-Fermi atom of the same radius becomes larger 
with Z . 

A difficulty arises when one seeks the next correction term to 
2 of Eq. (24) from Rijnierse's general solution for given by 
Eq. (Al) in the Appendix. In this case, Xj^ 2 is fixed by Eq. (25). 
However, a subseries in the square brackets of Eq. (Al) exists whose 
i-th term varies as e^x^+^ and consequently as x^ ^ at the isolated 
atom boundary. It is unfortunate that these terms contribute to the next 
correction to ^ because, in principle, one must sum an infinite series 
to obtain the corresponding coefficient. In fact, the first few terms of 
this subseries indicate at best a very slow rate of convergence. 

Jensen'* has determined the equilibrium atom radius of a spheri- 
cally symmetric and uniform distribution of electrons in the field of a 
nucleus by minimizing the Hamiltonian 

H = Jp^/^dT-e Jv^pdT + y / / "[r ^ r ' ^ dr , (27) 

where p is the electron density, is the nuclear potential and the con- 
stants are given by = (3/ 10) (S/tt) tt^ e^ a^ and ~ (3/4) (3/'^)^'^^ e^. 
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The first integral over the atom volume represents the kinetic 
energy, while the following two are the electron-nucleus and electron- 
electron potential energies, respectively, and the last integral 
represents the exchange energy. Mininizing the Hamiltonian with 

respect to the atom radius for the case of a uniform density P - Z/v 
in terms of the atom volume v = (4/3) i one obtains 




- 1-1 


1 + — ( 967 ^ 2 ) 1/3 22/3 


(28) 


which represents the atom radius determined by Jensen. Binomial 
expansion of this result yields agreement with Eq. (24) through terms 
of order 

In deriving Eq. (28), Jensen assumed that the electron distribution 
is uniform. He did not state any restriction on the atomic number Z. 
However, it is obvious from the argument of Gilvarry and that presented 
in I that Z->0 is the necessary and sufficient condition that the electron 
distribution in an isolated Thomas-Fermi-Dirae atom be uniform. It 
follows, in fact, that Eq, (28) is subject to the restriction that the 
atomic number Z be small, and that the result as expressed by Eq. (24) 
is no more accurate than the binomial expansion of Eq. (28) which 
corresponds to it. As a check, one notes that Eq. (28) implies that 
2 approaches a constant for large Z , which is clearly invalid since 
X|^ 2 is know to vary as Z^^^ in this limit,® 
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III. THERMODYNAMIC PARAMETERS 


For sufficiently small , the total Thomas-Fermi-Dirac energy 
with exchange, related to the pressure by = -du^j^/dv can be 
derived by integration of Eq. (7). Neglecting the resulting integration 
constant,® one obtains 

^ex=^o tl-2(l+32/3 e-i/5) (Po/pg)i/5]^ (29) 

in terms of the Thomas-Fermi energy Uq without exchange corresponding 
to Po > given by 


u„ = (6/5) ( 4 / 77 ) 2/3 Z'/Vx2 . (30) 

In analogy with the Thomas-Fermi-Dirac pressure p^^^ of Eq. (17), the 
energy u^^ can be expressed as a function of the Thomas-Fermi 
pressure and the corresponding energy Uj^ given by 




1 3^/^ (P„/p„)*/= 


(31) 


which one obtains by integration of Eq. (11). Accordingly, one has 


u = u„ 

ex M 




(32) 
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The last term in the square brackets of Eq. (3Z) is of form 
and has been retained since it is of the same order as ^ 3-t the 
isolated atom radius (as has been pointed out previously). The general 
solution for u^^ obtained by Rijnierse shows the presence of a term of 
form x^^ with coefficient independent of e , which clearly can be 
ignored. Hence, all terms have been retained in Eq. (32) for which 
are significant for = Xj^ ^ the case corresponding to the isolated 
atom, where Eq. (25) holds. Thus the ratio 

u^Vu = 1-2 (Pj/p)''/® (33) 

represents the asymptotic relation between the Thomas -Fermi-Dirac 
energy and the Thomas-Fermi energy u when the effects on the 
energies of not only the kinetic energy of the electrons but also the 
Coulomb interaction to first order are taken into account. This expres- 
sion should be valid for an extended range of low atomic number as 
compared to the corresponding result of I. 

The result of Eq. (33) should be valid in the two limits corresponding 
to infinite pressure for arbitrary atomic number and to vanishing atomic 
number for arbitrary pressure. For the limit p-*oo independent of Z, 
one obtains = 1; on the other hand, the last term on the right- 

hand side corresponding to the Coulomb correction vanishes for the 
limit Z-^0 independently of p, yielding the result of Gilvarry and as 
given in I. 
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If one defines the Thomas -Fermi-Dirac parameters 7 q and 
eo.ex exchange by v/ (To.ex " ^ ^ ~ " dlnPe/dlnv, 

respectively, in analogy with the corresponding definitions 
u= pv'/( 7jj-l) and = - dlnp/dlnv for the Thomas-Fermi ease 
without exchange, one can obtain and in terms of p^ for X|^ 

small. The results for and can then be expressed in terms 

of the Thomas-Fermi pressure Pj^ of March and the corresponding 
parameters y^^^ and determined by Uj, = Pmv/(7o,m“1) 

^0 M ~ “ ‘^fop,^/ dlnv, respectively. Consequently, one has 


^O.ex = ^0 




l-2(P„/p)‘/5-|z2/=(P,/p)2/5 


-1 

, (34a) 


^0,ex “ ^0 






-1 


(34b) 


where andp^^ have been replaced by their asymptotic 

equivalents 7 q, and p , respectively, for the case of small Xj^ . Mote 
that 7 q and in the present context are given by 


r„ =i. [i + (3'/Vio)xJ, 
3 


(35a) 


£„ =— [l+(3i/V20)xJ, 


(35b) 
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o 1 f) 

respectively, from the results of March, ’ and not simply by the limits 
7q = = 5/3 appropriate to the discussion of I. Thus, the parameters 

Jq and £(j are not identical for the case of low Z considered here. These 
expressions generalize the results of I to an extended range of low 
atomic number. 

The results obtained represent generalizations of the corresponding 
expressions first obtained by Gilvarry. For the case of vanishing atomic 
number independent of pressure p one cannot take P >> » in general, 

since is 0.01982Mb numerically. However, when the pressure is 
sufficiently high, a binomial expansion of Eq. (34) yields the approx- 
imations 




(36a) 


^0 , ex ^0 


(Pj/p)2/5 


(36b) 


generalizing the original results of Gilvarry. Since the Coulomb 
correction terms varying as in Eqs. (34) and (36) vanish 

for the limit Z-»0 at constant p, one notes that the results of I obtain 
in this case. 


The Thomas -Fermi-Dirac compressibility with exchange is 

defined by ySg = - dlnv /dp^^ in analogy with the Thomas-Fermi 
compressibility = - dlnv/dp without exchange. In terms of the 
differential parameters and introduced previously, one has 
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(37a) 




^^0 , ex ^ex^ 


-1 


^0 - ( ^0 P) ’ 


(37b) 


respectively, by definition. The compressibility ^o,ex (37a) 

when combined with Eqs. (17), (36b) and (37b) can be expressed as a 
function of Pj, and the corresponding compressibility m ~ ^ m Pm 
A ccordingly, for sufficiently small , one obtains 

-1 

, ■ (38) 


^0 , ex “ ^0 


1-4 (p„/p)‘"= 


upon replacing and (Sq by the corresponding p and /3 q , respectively. 
In this expression, one has 


^0 =|po~' tl + (3'/V5) x^] , 


(39) 


0 1 ft 

from results of March. ’ For the restricted situation where the 
pressure p is sufficiently high, the compressibility is approximated 
by 


'^O.ex " ^0 




(40) 
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Vhieh generalizes the result of I for the case of low atomic number. 
Note that, the compressibility reduces to the result of I for the limit 
Z -* 0 , independent of p • 

The results for thermodynamic parameters obtained in this section 
can be generalized on the basis of the results of Rijnierse in an obvious 
manner. 

IV, COMPARISON WITH NUMERICAL DATA 

The Thomas-Fermi-Dirac pressure as a function of the corre- 
sponding Thomas-Fermi pressure p has been plotted by Gilvarry for a 
number of values of Z from available numerical solutions for these 
quantities. The solid lines for constant Z in Fig. 1 where constructed 
from his results* The selection of values of Z includes the practical 
limits represented by Z = 92 and Z = 2. The limit line in the figure 
corresponds to the asymptotic expression of Eq. (19) for the limit 
Z - 0 . 

The asymptotic expression of the Thomas-Fermi-Dirac pressure 
Pex given by Eq. (17), which includes the Coulomb correction to first 
order, is shown by dashed lines in Fig. 1. One notes that at any value 
of p , the curves of constant Z approach the limit line as Z-» 0, in the 
family of solid curves corresponding to the exact solutions and in the 
family of dashed curves corresponding to the asymptotic result of Eq. 
(17). This fact is more obvious in Fig. 1 for the higher pressures than 
for the lower. Further, one notes that at sufficiently high pressure, all 


22 



curves (solid and dashed) approach the tangent line = p to the limit 
line. Thus, the march of the curves corresponding to the asymptotic 
expression of Eq. (17) is consistent with the limit line of Gilvarry and 
with the exact solutions. 

Figure 2 shows values of p^^ as a function of p of intermediate 
magnitude as solid curves on an enlarged scale, from numerical results 
of Metropolis and Reitz.^^ Corresponding curves from the asymptotic 
result of Eq. (17) appear dashed. Although the data are limited, it is 
apparent (most prominently for Z = 14) that the asymptotic expression 
approaches the corresponding exact result for fixed atomic number as 
the pressure increases. Further, it is clear from the figure that for 
a given value of Z , there exists a Thomas -Fermi pressure P , say 
P(Z)» such that for p > P(Z) the asymptotic expression P^^ of Eq. (17) 
is a closer approximation to the actual Thomas -Fermi-Dirac pressure 
than is the limit of Eq. (19) corresponding to Z-*0 • For example, 

P(14) is approximately 35 Mb while P(92) is about 350 Mb. It is 
obvious that for p > P(Z), Eq. (17) must be a valid approximation to 
the actual Thomas-Fermi-Dirac pressure for Z in the intermediate 
and large range rather than merely in the asymptotic limit Z-»0, 

The lower range of pressures is covered by Fig. 3, where the 
solid curves representing exact solutions for P^,,. vs P at low and 
fractional values of Z were plotted from results given by Gilvarry. 

The dashed curves yield the asymptotic results of Eq. (17). It is obvious 
that in the families of solid and of dashed curves, individually, the value 


23 



of the Thomas -Fermi-Dirac pressure approaches that given by the limit 
line as Z-»0 for fixed p or as p-*oo for fixed Z . Thus, the results are 
consistent in this sense with the theory given. Moreover, one notes 
that the courses of the dashed relative to the solid curves are essentially 
the same in Fig. 3 as in Fig. 2. Thus, one can infer that Eq. (17) should 
yield a valid approximation to the true Thomas -Fermi-Dirac pressure 
for low values of Z when one has p> P(Z). However, the data of Fig. 3 
are too limited to specify P( Z) numerically in the limit of low atomic 
number, although P ( 0.04 ) must be of the order of a megabar, as judged 
by extrapolation from the data. 

The salient feature of Fig, 3, however, is that for increasing values 
of Z at fixed p, the asymptotic values p^j^ of Eq, (17) become increasingly 
poorer approximations to the actual numerical values of p^^ at these 
lower pressures. Even for the lowest fractional value (0.04) of Z appear- 
ing, the value of p^^ corresponding to the limit line for Z = 0 is a better 
approximation to the true value than is that given by Eq, (17). This 
behavior obviously arises from the fact that the factor present 

in the correction term Z^^^ (Pi'^p square brackets of Eq, 

(17) causes it to overestimate considerably the Coulomb correction at 
low pressures where p < P(z)* a consequence, the truncated series 
in the square brackets of Eq. (17) yields a reasonable approximation at 
low pressure only in the mathematical cases of Z close to zero (for 
Z << 0.04 i at least). Hence, for the range of low pressures corre- 
sponding to Fig. 3, the asymptotic expression of Eq. (19) for the limit 
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Z — O yields a more useful approximation than does of Eq. (17) over 
the practical range Z > 1 starting with hydrogen. 

One notes that the asymptotic expression of Eq. (20) derived from 
the results of Rijnierse contains the further negative correction term 
Z4/3 j3/s within the square brackets, as compared to Eq. (17). 

Hence, the corresponding value of leads to a further overestimate 
of the true Thomas -Fermi-Dirac pressure. However, it can be shown 
that this trend is reversed by inclusion of the next correction term (not 
discussed in this work) from the results of Rijnierse, which would have 
the effect of adding a positive quantity to p^^ of Eq. (20). The necessity 
of such behavior was pointed out by Gilvarry in his initial discussion. 

It is clear that entirely similar statements can be made on the 
corresponding range of validity of the asymptotic expressions for the 
energy, the compressibility, and the thermodynamic parameters asso- 
ciated with the energy and the equation of state, as obtained. 

V. DISCUSSION AND CONCLUSIONS 

It is obvious that the extended asymptotic relation between the 
pressures on the Thomas-Fermi-Dirac and Thomas-Fermi atom models 
as obtained here is a direct consequence of the Thomas-Fermi-Dirac 
equation but has only a limited usefulness at low atomic number and 
low pressure. This circumstance is a result of the fact that the first 
(or first few) terms of an asymptotic series in Z for low atomic number 
in this case can yield a useful approximation for the limit of low atomic 


25 



number simultaneously with the limit of low pressure only if a fairly 
large number of terms is included from the results of Rijnierse. Note 
that this stricture does not apply to the expression obtained in I of this 
series for the limit of vanishing atomic number, because this result is 
exact (the atomic number Z does not appear in it). 

On the other hand, the asymptotic relation obtained between the 
pressures on the two models has an obvious use at higher pressures, 
since it specifies precisely how the pressure with exchange approaches 
the value of the limit line as the pressure increases for fixed atomic 
number. In this physicai case, a single correction term corresponding 
to the effect of the Coulomb energy reproduces the results from direct 
solution of the Thomas-Fermi-Dirac equation provided the pressure in 
question exceeds some limit dependent on atomic number, specified by 
P(Z) in the foregoing. 

The ultimate objective of this series of papers is to construct an 
asymptotic relation between the pressures on the Thomas-Fermi-Dirac 
and the Thomas -Fermi atom models which is valid for arbitrary atomic 
number and, to as great an extent possible, for arbitrary as well as 
high pressure. The results of this paper impose obvious constraints 
on the final solution of this problem. 
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APPENDIX 


In a recent paper, Rijnierse® has solved the Thomas-Fermi-Dirac 
equation by a method of successive integrations and shown that the 
solution \p can be obtained to any desired order for the case of small 
atom radius. Consequently, the Thomas-Fermi-Dirac pressure can be 
determined to any order, and Rijnierse has evaluated p^^^ to PqO ). 
The expression of Rijnierse for the pressure is 

Pex “ Po -3^^^ x^/4 - 5eXj^/4-3^'^^ - Rj x^ 

+ ^2 (4-3^'^^/xj^) + R 3 ex^] (Al) 

in terms of the constants 

_ 54,929 _ _ 31,500 ^ ./g 38,880 

^ 907,200 ’ 2 ^ 907,200 ’ ^ ' 907,200 ’ (^2a, b, c) 

In analogy with the Thomas-Fermi-Dirac pressure of Eq. (16), the 
pressure of Eq. (Al) can be expressed as a function of the corre- 
sponding Thomas -Fermi pressure Pg given by 

pR = Po x^/4 - Ri x^ + R 2 x^ In (4-32/Vxj^)] , (A3) 
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determined from Eq. (Al) by setting e to zero. By making Pg a factor 
of Eq. (Al) and adding and subtracting appropriate terms, one obtains 

p^^ = Pjj [l-(5ex/4-3i/3) (1+3^/3x^20 + 35/3x2/400) 

- (6x^4) (1 + 31/3x^10) - (7-31/3/160 -R 3 ) ex3] , (A4) 


retaining only terms of significant order of approximation. This 
expression enables one to write p^^^ as a function of Pg by 

p.. = Pr [1-(Po/Pr)‘'®- Z"''* (P/Pr)"'= - Z*'’ (A5) 

in terms of the constant pressure Pj of Eq. (21) in the text and the 
constant pressures and P^ previously defined. 
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FIGURE CAPTIONS 


Fig. 1. Curves of the pressure with exchange as a function of the 
Thomas -Fermi pressure p for selected values of atomic number Z. 
The limit line corresponds to Eq. (19) while the remaining solid 
curves correspond to numerical solutions. For comparison, the 
dashed curves correspond to the asymptotic expression of Eq. (17). 

Fig. 2. Curves of the pressure with exchange as a function of the 
Thomas -Fermi pressure p for selected values of atomic number Z 
The limit line corresponds to Eq, (19) while the remaining solid 
curves correspond to numerical solutions. For comparison, the 
dashed curves correspond to the asymptotic expression of Eq. (17). 
The scales correspond to the intermediate decades of p^^ and p in 
Fig. 1. 

Fig. 3. Curves of the pressure p^^^ with exchange as a function of the 
Thomas-Fermi pressure p for selected values of atomic number Z 
The limit line corresponds to Eq. (19) while the remaining solid 
curves correspond to numerical solutions. For comparison, the 
dashed curves correspond to the asymptotic expression of Eq, (17). 
The scales correspond to the lowest decades of p^^ and p in Fig, 1. 
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